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ABSTRACT 

A family of simple, periodic basis functions with “built-in” discontinuities are introduced, 
and their properties are analyzed and discussed. Some of their potential usefulness is illus- 
trated in conjunction with the Fourier series representation of functions with discontinuities. 
In particular, it is demonstrated how they can be used to construct a sequence of approxima- 
tions which converges exponentially in the maximum norm to a piece-wise smooth function. 
The theory is illustrated with several examples and the results are discussed in the context 
of other sequences of functions which can be used to approximate discontinuous functions. 
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1. Introduction. Fourier series are used widely in many branches of applied 
mathematics. For example, they are often used together with separation of variables 
to construct analytical solutions to boundary value problems for differential equations, 
and with a variety of spectral methods to find approximate solutions to these prob- 
lems numerically. For practical purposes, approximate solutions to these problems are 
often obtained using only a finite number of the terms in a Fourier series. This trim- 
cation procedure may lead to nonuniformly valid approximations. In particular, when 
the function being approximated has a point of discontinuity, the Gibbs phenomena is 
present. It is well known (see, e.g., [1] or [12]) that the magnitude of the “overshoot” 
associated with the Gibbs phenomena is not eliminated by increasing the number of 
terms in the approximation, but rather the overshoot approaches an asymptotic limit 
of about 18% as the number of terms approaches infinity. In addition, for any Fourier 
series partial sum, the “oscillations” caused by this phenomena typically propagate into 
regions away from the singularity, and, hence, degrade the quality of the partial sum 
approximation in these regions. 

In a series of papers, Gottlieb and several of his co-investigators [3], [4], [5], [6], [7] 
have proposed a way of overcoming the Gibbs phenomena. Their technique involves 
the construction of a new series using the Gegenbauer polynomials. For a fimction 
/ that is analytic on the interval [—1,1], but is not periodic, they prove that their 
technique leads to a series which converges exponentially to / in the maximum norm. 
To do this, they require that the parameter A, which appears in the weight factor 
(1 — rc 2 ) A-ly/2 , grows with the number of Fourier modes considered. Recently, Geer 
[2] introduced and studied a class of approximations {F NM } to a periodic function 
/ which uses the ideas of Pade, or rational function, approximations based on the 
Fourier series representation of /, rather than on the Taylor series representation of /. 
Each approximation F'n,m is the quotient of a trigonometric polynomial of degree N 
and a trigonometric polynomial of degree M. The coefficients in these polynomials are 
determined by requiring that an appropriate number of the Fourier coefficients of Fm,m 
agree with those of /. Explicit expressions were derived for these coefficients in terms of 
the Fourier coefficients of / and it was proven that these “Fourier-Pade” approximations 
converge point-wise to (/(a: 4- ) + / (x~))/2 more rapidly (in some cases by a factor of 
1 jk 2M ) than the Fourier series partial sums on which they were based. Although these 
approximations do not “eliminate” the Gibbs phenomena, they do mitigate its effect. 
In particular, the asymptotic value of the magnitude of the overshoot is reduced to 
about 6%, and, outside a “small” neighborhood of a point of discontinuity of /, the 
“unwanted” oscillations can (for practical purposes) essentially be eliminated. 

To fix notation, we let / be a 27r-periodic function with enough smoothness and 
regularity properties so that its Fourier series exists and converges to ( f(x + ) + f(x~))/2 
for all —7 r < x <7 r. Let C M [— 7 r, 7r] denote the class of 27r-periodic functions that have 
at least M continuous derivatives on [— 7r, tt[. Then it is well known (see, e.g., [1]) that, 
if / e C M [— 7r, 7r], then the Fourier coefficients {a.j,bj} of / are 0(1/ j M+2 ) as j — * oo. 
Moreover, if / is 2n— periodic and analytic on [ — 7T, 7r] , then there exists a constant 6, 
with 0 < 6 < 1, such that aj and bj are O(0 J ), as j — ► oo. In this case, it follows that 
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the Fourier series of / converges exponentially to / in the maximum norm. 

In this paper, we introduce a new, simple class of periodic “basis” functions which 
have certain “built-in” singularities, and which can be used to construct a sequence of 
approximations which converges exponentially to / in the maximum norm. In particular, 
this implies that the Gibbs phenomena can be completely eliminated, even when / has 
several points of discontinuity in the interval [— tt, 7r] . In section 2, the first, two o 
these basis functions are introduced and it is shown how they can be used to improve 
the convergence properties of a function with one or more singularities in the interval 
f-7r,7rl. These ideas are illustrated by two simple examples in section 3. In section 4, 
the complete class of basis functions are introduced, and various convergence properties 
of approximations constructed using them are studied in section 5. In particular, it 
is proven that, under only mild restrictions on the original -function /, they can be 
used to construct a sequence of approximations that converge exponentially to / m 
the maximum norm. The results of the theorems of section 5 are illustrated by three 
more examples in section 6. We discuss our results in the final section and compare our 
approximations with those used by other investigators. 

2. Basis functions S 0 and Si. Consider the 27r-periodic functions 


( 1 ) 


So(x) = 


sin(x) 


2\/2 y/l - cos(x) 


and Si(x) = — cos(x), 


which are illustrated in Figure 1. Using the Taylor series expansions of sin(x) and cos(x) 
about x = 0, it is straightforward to show that 


(2) 5„(x) = ^{^+0(x 2 )}, S 1 (x) = w{i+0(x J )}. as x — * 0. 

From these relations it follows that S 0 and S 1 satisfy the following jump conditions at 


x — 0 : 

( 3 ) 


[So(0)j = 1, 


dS o 
dx 


( 0 ) 


= 0, [Si(0)]=0, 


'dSi 

dx 


( 0 ) 


= 1 . 


Here we have used the notation [<?(*)] = g(x + ) - g(x-). (In the second and fourth 
equations in (3), the derivatives of S, axe understood to be derivatives from the right 
and “from the left” , respectively.) We also note that both So and Si have continuous 

derivatives of all orders for -tt < x < 0 and 0 < x < n. 

Consider now the problem of approximating the piece-wise continuous (smooth) 

27r-periodic function f{x), which can be represented at points of continuity of / by its 
Fourier series, i.e., 


N 


f(x)= lim /*">(*), F^H^) = T + S a J cos 0' :c ) + ^ sin ^ ) ’ 

J x 7 N—> oo ^ j— 1 



= 1 f m 

7T J — 7T 



^ dx, j = 0,1,... 


( 4 ) 


We shall assume that f(x) and/or f(x) may have discontinuities at a finite number 
of points, say at x = Xj, j = 1,2, ...,n, in the interval (—■ 7r,7r], and that the one-sided 
limits of / and /' exist and are finite at each of these points. 

To construct a more rapidly converging Fourier series type representation of /, we 
first define the function 

71 

(5) Ri(x) = {4), AO* “ x i) + a i,j S i{ x ~ x i)} » 

i = i 

where the constants {A 0j -} and {Aij} are determined so that the “jumps” in R\(x) 
and its derivative coincide with the corresponding jumps in f(x). Using the relations 
(3), this requirement leads to the expressions 

( 6 ) A o,j ~ [f( x j )] . Aij = [f(xj )\ , j = 1 , 2, ...n. 


Then we write 

(7) f{x) = Ri{x) + fi(x), where f 1 (x) = f(x)-R 1 (x). 

By the way it has been defined, the function fi(x) is C 1 [—tt, 7t] (at least) and, hence, 
its Fourier series will converge at a faster rate than the Fourier series of /. To construct 
its Fourier series, we first note that Sq and S\ can be expanded in a Fourier series as 


( 8 ) 

and 

( 9 ) 


S 0 (x) - 


1 sin(x) ^ 4 j 

= 2^ b Qj sm(jx), boj = 


2v2 J 1 — cos(x) j= i 


7 r 4jf 2 — 1 ’ 


<^1 o 00 

5i(x) = \/l -cos(rr) = + £ a U cos 0' x )> a U = 


j = l 


tt( 1 - 4 j 2 )' 


Then, denoting the Fourier coefficients of /i by {a^,6^}, we find, using equations 

(4)-(9), 

TL 

a'p = a.j - ]T) {A 1>s a hj cos (jx s ) - A 0)S b 0 ,j sin(jx s )} , 

S=1 


(10) bf ] = bj - An(jx s ) + A 0tS b 0J cos(jx s )} , j = 1,2,.... 

5=1 

As we shall demonstrate explicitly in the next sections, the coefficients {a^, b^} decay 
faster (by at least a factor of 1/j 2 ) as j — » oo than the original Fourier coefficients 
{a.j,bj}. Using these definitions, we define the new approximations 

(i) N 

(11) x ) = R\(x) + -2- + ^2 cos (jx) + b^ sin (jx), N > 0. 

2 j = i 
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These functions are the first in a sequence of “better” approximations to /, which we 
shall illustrate, generalize, and then analyze and discuss in the following sections. (We 
note that the general idea of “subtracting” from the Fourier series of / the Fourier series 
of another “simple” function with an appropriate singularity has been considered by 
other investigators, for example, by Tolstov [13], pp. 144-147. However, this idea does 
not seem to have been pursued in any systematic or exhaustive manner.) 

3. Examples. Before we generalize the family of approximations } and 

discuss some of their convergence properties, we illustrate these approximations by two 
examples. 

Example 1: We consider first the function f(x) defined by f(x) = x + 7 r, for 
— tt < x < o, f(x) = x — it, for 0 < x < 7r, and f(x + 2-n) = f{x), for all x. Then 
x = 0 is the only point of discontinuity of f(x) in the interval -tt < x < tt. Thus, we 
set n = 1 with xi = 0 in equation (5) above and use the facts that [/(0)] = —2n and 
[/'(0)] = 0 in equations (6) to define 

7r sin(x) 

(12) tfiW - 4uS>(*) + = -y| ^_ cai{x) - 

Thus, Ri(x) has the same jump at x = 0 as the function f(x), which is clearly evident 
in Figure 2, where we have plotted both f(x) and Ri{x). Hence, we now consider the 
function /i(x) = /(*) - Ri(x), which is C 1 ]- tt.tt], and is also shown in Figure 2. Using 
equations (10), we expand f\(x) in a Fourier series as 

00 2 

(13) /i(*) = £ 6 ? )sin 0s)> b i = j j 2 _ 1) • 

We note that the coefficients {bf } } are 0(l/j 3 ) as j — * oo, while the Fourier coefficients 
b = —2/j of / are only 0(1/ j) as j -+ oo. Then, using equation (13) in equation (11), 

we define the approximations 

(14) f {l ' N) (x) = Ri{x) + J^bf ] sm(jx). 

In Figure 3 we have plotted both f(x ) and / (1 ’ 3) (x), where the good agreement is 
evident. 

Example 2: As a somewhat more complicated example, we consider the function 
/ defined by 

— (x + 7r/3), — tt < x < — 7t/3, 

2x, — -tt/3 < x < 0, 

1, 0 < x < 7t/2, 

(x — 7r/2 ) 2 + 1/2, tt/2 < x <i r, 

and f(x + 2tt) = f(x) for all x. This function is illustrated in Figure 4, along with the 
Fourier series partial sums F <N) {x) of f(x), defined by equations (4), for N = 10, and 
the error (difference) f(x) — F^ 10) (x). 
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The function / has discontinuities at the n — 4 points x x = — 7t/3, x 2 = 0, x 3 = tt/ 2, 
and x 4 = 7T. Using the definition (15) and the relations (6), we find that the coefficients 
{^4 0 j} and {^4i,j} in the definition (5) of Ri(x) are given in Table I. 


Table I: Coefficients Akj for Example 1. 


3 

A oj 

Aij 

1 

-2tt/3 

3 

2 

1 

-2 

3 

-1/2 

0 

4 

(8tt — 6 — 3tt 2 )/12 

-(1 + tt) 


The functions /(x), Ri(x), as well as the difference fi(x) = f(x) — Ri(x), are shown 
in Figure 5. Obviously, /j(x) is much smoother than /(x), and in fact, the three-term 
Fourier series partial sum of fi(x) is virtually “identical” to /i(x), to within plotting 
accuracy. Thus, the approximation /( 1,3 )(x), which is shown in Figure 6 along with /(x), 
is a much better approximation to /(x) than the original Fourier series representation 
of /, even when many more terms are included in F^(x) (cp. Figure 4). 

4. Generalizations. The basis functions introduced in the section 2 can be gen- 
eralized and extended to allow us to construct functions which have appropriate dis- 
continuities at, say, the n th derivative. To this end, we define the functions S n (x) by 

nk—3/2 

s 2 k (*) = sin(x) (1 - cos(x))*- 1/2 , 

nk-l/2 

(16) (l-cos(x)) t+1/2 , *= = 0,1.2,.... 

Then it is straightforward to show that £>„(x) is C n ~ l [— tt, 7 r], while the jump in its n th 
derivative at x = 0 is 1. More generally, each S n has continuous derivatives of all orders 
for — 7 r < x < 0 and 0 < x < 7 r, and, at x = 0, satisfies the jump conditions 

[S<’>(0)] s = 0, if 9 < n; [S<">(0)] = 1; 

L\ ax / J x=o 

[s! 2 ' +1 >(0)] = [sj 2rt (0)] = (-1/4)', p > 0; 

[-sg'A(o)] = [4' +1> (o)] = 0 , p > k > 1 , 

(17) 





In particular, specific values for the jumps of [sg p) {0)} are summarized in Table II for 
a few small values of k and p. 


Table II: The jumps = [-YjY (°)j 



p = 0 

1 

2 

3 

4 

5 

6 

k = 0 

1 

-1/4 

1/16 

-1/64 

1/256 

-1/1024 

1/4096 

1 

— 

0 

1 

—5/2 

91/16 

-205/16 

7381/256 

-33215/512 

2 

0 

0 

1 

-35/4 

483/8 

-12485/32 

631631/256 

3 

0 

0 

0 

1 

-21 

2541/8 

-68497/16 

4 

0 

0 

0 

0 

1 

— 165/4 

18447/16 

5 

0 

0. 

0 

0 

0 

1 

-143/2 

6 

0 

r 0 

0 

0 

0 

0 

1 


Using the definitions (16), we find that the Fourier series of S n can be expressed as 


oo 00 

(18) S 2 fc+i(z) = + 12 a 2k+ijCos{jx), S 2k {x) = sin 0' x )> 

1 7 = 1 j = 1 


a 2k+l,j 


= (-l) 


k - hi 


4 fc+1 


7 r 


n (4j 2 — (2i + l) 2 ) 


i = 0 


(19) 

for k = 0, 1. 2, ... 


62 kj — — j&2k+l,j > 3 2, ... , 

(For convenience in certain formulas which follow, we also define 


^2k,j 


= b 


0, for k > 0.) 


— l/2A:+1 } j — v/ , ivi • * ^ ■ j ,,, 

Now let the 27T-periodic function / have possible discontinuities in / for 0 < k < 
M at x = x , j = 1,2, n, where -tt < x 3 < 1 r. (Here we assume that all of the 

one-sided limits f (k) ( x j ) and f {k) ( x j) exist and are finite -) We define the functlon 


M n 

(20) Rm { x ) = F AkjSkjx — Xj), 

k=0j=l 

where the constants {A fcj } are determined so that /m(x) = f(x)-R M (x) is C M [— 7r,7r]. 
Using the relations (17), we find that these constants can be determined recursively from 

the relations 


( 21 ) 


z— 0 


Thus, the function / M (x) will be C M [-7r,7r], at least, and hence its Fourier series will 

converge at a faster rate than the Fourier series of f . 
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Once Rm(x) has been defined, we define the family of approximations to / 

by 

AM) N 

f( M ’ N )(x) = R m (x) + -9— + jo af q cos {jx) + bf 1] sin(jx), 

z j=i 

M n 

aJ M) = A k>s {a kj cos (jx s ) - b kJ sin {jx,)} , 

k=zO 5=1 

M n 

(22) 6j M) = bj-^2Yl Ak ’ s ( a *d sin U x *) + cosO'x s )} , j = 0, 1, 2, ... . 

fc=0 s = 1 

These approximations have several nice convergence properties, which we shall discuss 
in the next section. 

Returning to Example 1, since the function fi(x) (see equation (7)) is C 1 [— 7r,7r] 
and satisfies the jump conditions [/"( 0)] = -n/2 and f/"'(0)] = 0, we use the definition 
of Rm { x ) i n equation (20) with M = 3 to define 

i? 3 (x) = -2irS 0 (x) - S 2 (x ) 


7 r sin(x) 7r r — 

(23) =-^7T^w-^ sm(l) ' /r ^ cs(l) - 

Then / 3 (x) = f(x) - R 3 (x) is C 3 [- 7 r, 7 r], and the Fourier series coefficients in the 
definition (22) of / (3>iV) (x) are given by 

18 

( 24 ) aj } = 0, j > 0; b ( j } = ~ j{4]2 _ i)(4j2_ 9 )’ 3 - h 

We note that bf ] = 0(1/ f) as j -* oo, while the original Fourier coefficients of / are 
only 0(l/j) as j — ► oo. 

5. Convergence properties. We now investigate some of the convergence prop- 
erties of the function defined in the previous section, as M, N — * oo. For this 

purpose, we first define the errors 

E^ m ' n \x) = /(x) - / (M>A °(z) 

OO 

(25) = Y/ cos(jx) + sin(jx), M > 0, N > 0, 

j=AT+l 

where the coefficients {a { ^\bf I] } are defined in equations (22). Since the function 
f M (x) = /(x) - Rm(x) is C M {- 7T, tt], it follows that both af n and b\ M) are 0(j“ (A/+2) ), 



as j -» oo (see, e.g., [1]). Thus, there exists a constant K = K{M,N ), independent of 
j, such that 


(26) 


a 


(M) 


< 


K 


\b™\ < 


j I — jM+ 2’ 


K 

am- 1 - 2 ’ 


for all j > N + 1 . 


Using the bounds (26), we find from equation (25) that 


(27) 


oc 

£ (m ’ n) (x)| < 2K 53 


1 


„'M+2 


< 


2/^ 1 
M+17V M+1 ’ 


which follows easily from the integral comparison test. 

Suppose now that K satisfies a bound of the form 

(28) K < K 8 m M p M\ , as M, iV — » oo, 

where K and 6 are constants independent of M and N. Then, if we let N be propor- 
tional to M, say, N = AM, where A > d/e, it follows that the errors E [M ' N) {x) decay 
exponentially to zero as M — ► oo. To see this, we substitute (28) into (27) and use 
Stirling’s approximation to the factorial to write 

2 K 1 ^ - M p 9 M V2^MjM/ef_ 

M + 1 iV M + J - K M + 1 (A M) M+l 

(29) < . 

where K is another constant, independent of M and N, and (9/ Ae) < 1. Thus, we have 
proved the following theorem. 

Theorem 1. Let f(x) be a piece-wise smooth, 2n -periodic function with possible 
discontinuities in f and/or some of its derivatives at a finite number of points in the 
interval — n < x < tt. Suppose that the Fourier coefficients {a) , b } } defined by 

equations (21 )-(22) satisfy the inequalities (26), where K satisfies the bound (28). Then 
the approximations {/(".*")} defined by equations (20)- (22) converge exponentially to 
f(x) in the Loo— norm as M — > oo, if A is chosen so that A > 9/e, i.e. 

max \f(x) - /<"'“>(*)! < KW-W (/ Y , 

-TT<x<Tr J 



where K is a constant independent of M. 

To illustrate this theorem, the coefficients {^ Mj } for Example 1 for M = 2m + 1 
are easily shown to be given by 


(30) 


^(2m-hl) 


(1 3-5---(2m+l)) 2 
3 n”oW 2 -(2i+l) 2 )’ 


m = 0, 1, 2, ... . 
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Noting that 1 • 3 • 5 ■ • • (2 m + 1) = (2 m + l)!/(2 TO m!), we use Stirling’s approximation 
to the factorial in equation (30) to write 


( 31 ) 


, (2m+l) 
°3 


< K 


r,VMM\- eU 


,M+2 


9 m = 


2 M n^o(l-((^ + 1 )/2j) 2 )’ 


where K is a constant independent of j and M. In the Appendix we show that there is 
another positive constant, K , such that 

(32) e M < ke 0 M , e 0 = e 0 (\) = 1/ (2^/1 - i/(4A 2 )) , 


for all A > A, where A is any arbitrary number greater than 1/2, and for all j > 
N + 1 = AM + 1. Thus we see that bf m+V) satisfies the bound (26), where K has the 
form of the right side of equation (28) with p — 1/2 and 6 = $o(A). (See the Appendix, 
also, for a derivation a sharper bound for 6.) Then the condition that A > 0 o /e leads 
to the requirement that A > A 0 = \/l + e 2 /2e = 0.533. In Figure 7 we have plotted 
the normalized errors M (Ae j 6q) m (x) as a function of x, with A = 1 (and hence 
N = M) for M = 5, 9, and 13. We note that the error, when regarded as a function of x, 
is much more uniformly distributed over the interval [— 7r,7r] than is the corresponding 
error for the original Fourier series, which is highly concentrated around x = 0. In 
Figure 8 we have plotted the normalized L ^— errors M(\e/0 o ) M max |£^ M>jV) (x)|, with 
N = AM as a function of 1/M for several choices of A > Ao- The figure clearly illustrates 
that these quantities are bounded as 1/M -> 0 and, hence, that the errors E (M ’ N) {x) 
decay uniformly and exponentially to zero as M increases. 

To investigate the convergence properties of the approximations for a 

general function /, we first consider the following lemma. 

Lemma 1: Let f(x) be a 2 tt— periodic function, n a positive integer, and the points 
X 0 = -7T < Xi < ... < x n — 7T be defined such that fix) = Pk{x) for x fc _! < x < x k) k = 

1,2, ..., n, where p k is a polynomial of degree at most m. Then the Fourier coefficients 

of f defined by equations (4) can be expressed as 

1 n f /[ m /2] /'i'ifc+1 r A 

o 3 = - E | sin 0' x ^) ( E - ? - 2 fc+ r / {2fc) (*-)] ) 

71 s=l ( \ fc=0 J } 


+ cos(jx s ) 


1 n 
bj = - V 


+ sin (jx s ) 




where, m the upper limits of the inner summations, we have let [q] denote the greatest 
integer not exceeding q. 

The proof of this lemma follows easily by using the formulas (4) and integrating 
each of these expressions by parts m times (see, e.g., [10], pp. 489-493, for more details). 

Inserting these expressions into equations (4) and rearranging the resulting series, 
we can express the Fourier series of / as 

1 n n 771 

- (f(x + ) + fi x ~j) = T + [/ <fc) ( a; *)] h k( x ~ X *)’ 

2 v z s=lfc=0 


(33) 


h 2 k{x) = 


(-1)* £ 


sin (jx) 

j2k+l ’ 


h 2 k+l( x ) 


(-l)* +1 g 


cos (jx) 

~j2h+2 ’ 


for k = 0, 1, ... . Hence, to establish the exponential convergence of our approximations 
for a “piece- wise polynomial” function /, as defined in Lemma 1, it is sufficient to 
establish that our approximations converge exponentially to each of the functions hk(x) 
defined in equations (33). In fact, we shall show that 


(34) 


, , x h {M,\M)( s I < Ko M* / 0o\ 
ft,(x) - h\ (*)| < y-r ^ eA ) . 


o.i,... , 


and hence 
(35) 


max 

-7T<X<7T 


hjx) — h 




(x) 


< 


Konf 
M i\ 



i = 0,1,... , 


where K 0 is a constant, independent of M, 9 0 = 1/ (2^1 - 1/(4A 2 )) , with A = N/M , 
and 


(36) 


M N 

h% M (x) = E A^S^ + ^bf^smUx), 

k=2p j=l 

k even 


(37) 


M 

«f(x)= Y, 4 2p+1, s*W + 


a 0 


fc=2p+l 
k odd 


N 

+£4 

i=i 


(2p+l,M) 


cos(jx), 


for p = 0,1,... . Here the constants 4 2p) 3nd 4 2p+ °> 35 wel1 35 the coefficients {&f p,M) } 

and {a^ 2p+1,M ^}, are determined as described in section 4. 

To^show that equations (34)-(37) hold, we shall use induction. The function h Q = 
(— l/ 27 r)/(ar), where f(x) is the function of Example 1. Hence, using equation (33) with 
k = 0 and equation (36) with p = 0, we can write 


h 0 (x) - h 


( M,N ) 
0 



e 

j—N+l 


sin (jx), 


10 


(38) 



where bj M \ with M — 2m + 1, is given by the right side of equation (30). Using the 
bounds (26) and (28), with 6 = 6 Q , p = 1/2, and N — AM, as well as the remarks below 
equation (30) and Stirling’s approximation to the factorial, we can write 


h 0 (x) 


h (M,N) 

n 0 


(*) 


< 


K o 
M 



J 


where K 0 is another constant, independent of M. Hence, equation (34) holds with i = 0 
and the exponential convergence of our approximations to ho is established. 

Suppose now that equations (35)-(37) hold for i = 0, 1, ...,q — 1. To investigate the 
behavior of h^ M,N \ we observe first that 

(39) K+i( x ) ~ h i( x ) and S 2 i+i( x ) = S 2i (x), 

for all nonnegative integers i. (Here the primes denote differentiation with respect to 
x.) Then, using the first relation in (39), with i = q — 1, we can write 

h q (x) = h q (0) + f h' q (t)dt = h q { 0) + f h q ^(t)dt 

J 0 J 0 


(40) = h[ M ’ N) {x) + E^ N \x), 

where 

hf<^(x) = h,(0) + f hf_[ N \t)dt, 

J 0 


(41) E[ m ’ n \x) = r {Vi(t) - h^i N \t)} dt. 

J o 

Using our induction hypothesis and equation (34) with i — q — 1, we find from equations 
(40) and (41) that 


h,(x) - h^(x) I = I £<»■“> (x)| < | ViW - * 

J 0 


and hence 
(42) 


*1 


-i /a.\M Ko\x\U6o\ M 


*\Kq t*- 1 (e^ 

M(q- 1)! leA 


dt < 


M q\ \e\ ' 


- h^(x ) | < (^) 


M 


Thus, it just remains to be shown that h ( q M ' N> (x), defined in equations (41), has the 
form indicated in equations (36) and (37). To see that this is the case, if q is odd , we 
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use equation (41), equation (36) with 2p = q-l, and the last relation in equation (39) 
to write 


M JV 

h( MN \x ) = >>,(0) + f \ At i) s k (t)+Y,t < r' M)sin w 

q J0 k=q - 1 3 = 1 


, k even 


> dt 


N 


Af 


= M 0) + E(i/j)^" 1,m) + E 4 9 ' 1) ^ +1 (x) + E(-i/j>r ' )cos (^)> 


j=l 


k=q— 1 
Jc et?en 


N 

E 

j=l 


which has the form of equation (37) with 2 p +1 q , where 


JV 


a^ M) =2^ 5 (0) + ^(l/j>)' 

j=i 






(43) 4 9) = 4-4 a t M) = i-Vm 

If q is even, we use equation (41), equation (37) with 2p + 1 = q - 1, and the fact 
that h q ( 0) = 0, to write 


M 


(g-l,M) N 

h w.m {x) = £ J £A ( r l) S t (t) + +E4 , " , ' M) cos( 3 'i) 

k odd 


> dt 


N 


(44) 


= ^(x) + E(l/j>5 , - 1 ’ M) sin0rr), 


J = 1 


(45) 


»,(*) s f; 4 ,_1) f &(*)* + 

i i d Q 


c4 9_1 ’ M) X 


k—q— 1 
fc odd 


We now make some observations concerning the function First, we note that h q is 
an odd function of x, with h q ( n) = h q (— - k) = 0. The fact that h q is odd follows directlj 
from equation (45) and the fact that each S k is an even function of x whenever the 
index k is odd. To see that h q vanishes at x = ft, we note first that, by letting A — * oo 
in equation (37) with 2p + 1 = q - 1, we can write 


jVf (qf— l,Af) oc 

v iM= E +E4'' Mn “ s « i )- 

fe=7-l J =1 

odd 
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Integrating this expression from x = 0 to x = n and using the definition of 
(equation (33)), we find 

(9-1, AT) m „ 

(46) = - E A k / S k {x)dx. 

Z k=q - 1 ' /0 

k odd 

Setting x — ir in equation (45) and then using equation (46), we find that h,_i(7r) = 0. 
Using this first observation, we note next that h q -\ has q — 1 continuous derivatives for 
— 7 r < x < 7 T, and that the only discontinuities in higher (even) order derivatives occur 
at x — 0. Thus, using the ideas of section 4, we can write 

M OC 

(47) V i(z) = ^ i*,5 fc (a:) + ^6j M) sin(ja;), 

k=q j—l 

k even 

where the constants Ak are determined recursively by 

Aj = Y. A k~ l) [ 5 it J_1) (0)] - [4 J) (0)] > even, 

k~q—\ k=q 

k odd k even 

and the coefficients b Ar} are determined by the right side of the last equation in (22), 
with x s — 0, s = Ak, and 


E 4 

k=q - 1 
k odd 




k=q— 1 
k odd 


Thus, when the right side of equation (47), with the upper limit of the second summation 
replaced by N, is substituted for h q --y(x) in equation (44), we see that h ( q M ' N) has exactly 
the form of the right side of equation (36) with 2p = q. 

We now summarize these results in the following theorem. 

Theorem 2. Let f be a 2ix -periodic, "piece- wise polynomial” function, as defined 
in Lemma 1 . Then the sequence of approximations f( M AM) ^e/med in section 4 converge 
eocponentially to /, with 

4(sf' 


where 

n m 

= + £ £ [/“M ^"'"V - *.) 

Z 5=1 fc=0 

and K is a constant (which depends on the particular function /, but is independent of 
M ) such that 
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where K 0 is the constant which appears in the bounds (35). 

Finally, we consider the case when / (x) is a 27r-periodic, piece-wise analytic func- 
tion. We let the n points x 0 = -ir < x x < ... < x n = tt be defined such that 
f( x ) = ft( x ) for Xi -! <-X<Xi,i = 1,2, ...,n, where each /,(x) is analytic on the closed 
interval x,_i < x < x„ and such that there exists a point on, with x*_i <oq< x it such 

that (x, Xi- 1) < Pi, where p, is the distance between ex, and the nearest singularity of 

ft in the real or complex plane. Then, by Taylor’s theorem, we can write 


/i(x) — Pi,m{x) + , 


(48) Pi,m(x) = 53 


fPifli) 


°° f^(a) 

(x-a+, Ri,m( x )= 53 -~j (x-ai) 3 . 

j=m+\ J ' 


3=0 


We now assume (see, e.g., [9]) that 


(49) 


/i 0 , (Oi)| < C,jj< l<i<" 


where Ci is a constant, independent of j and />, Using equations (48) and (49), we can 
write 


max | fi(x) — Pi,m{ x )\ — l-^vnO 1 )! — Ci 


( 50 ) uu ~ / 


Xi—\ 

k Pi 


m+1 


where C t = C lPl /{ Pl - Xi + a,) and (x t - x i _ 1 )/p l < 1, from our definition of the points 
{x,} . We now define 

- ( X, Xi — | 

(51) 


C = max Ci, p = max 

l<i<n !<*<” 


< 1, 


SO that from equations (50) and (51) we can write 
(52) 


max |/t(x) — Pi,m( x ) I ^ Cp 

X t _l <X<Xi 


, 771+1 


We now define g{x) as the 27r-periodic, piece-wise polynomial function with g(x) - 
P im (x) for Xi-i < x < x^ i = 1,2, ...,n, and also define the approximation 




(x) = ^+tf:[9 {k) ^)] h { k M ’ N \ x -x s ) 


5=1 k = 0 


where h[ M ’ N) (x) is defined in equations (36)-(37). Using Theorem 2 and equations (33) 
and (35), we see that 


| S (x) - </"■«>(*) 


< Y, Y |[s W (*>)] (M* - *<) - hf- XM) (x - x .)) 


14 



( 53 ) 





We now assume that the jumps in f (k) satisfy a bound of the form 


(54) 


f {k) (x s ) 1 <A£ k k\, 5 = 1,2, 


where A and f are positive constants, independent of k and s. (Here we may assume, 
without loss of generality, that fyr > 1.) Since the jumps in g {k) are arbitrarily close to 
the corresponding jumps in f^ k \ we assume that equation (54) also holds for the jumps 
in g^ k \ and then use this bound in equation (53) to write 


(55) 


max 

7r<x<7r 


9 (x) x) 


<nA% 

M 





We now combine the bounds 


max 

— 7T<X<7T 


| /(*) 




(x) 


< 


(52) and (55) to write 

max \f{x)-g{x)\+ max | g{x) - g {M ' XM] {x) 

-7r<£<7T -7T<X<7T I 


(56) 


< Cg 


m+1 + nA^ 

M 


'(tt £) m+1 - 1 
7r£ — 1 



We now set m — M in equation (56) and note that the second term in equation (56) 
is dominated, for large values of M, by a constant times where ( = rr^Bo/eX. 
Since we want the magnitude of this quantity to be less than 1, we must require that 
A > fye 2 + 47 r‘^/ 2 e. Then, with this stipulation, we find from equation (56) that 


(57) 


max 

-7T<£<7T 


/(x) - <x)| < k 


as M 


oo. 


where K\ and K 2 are constants, independent of M, while 0 < g < 1 and 0 < £ = 
7r(0 o /eX < 1. Thus, we have proved the following theorem. 

Theorem 3. Let f be a 2 tt - periodic, ‘‘piece-wise analytic’ function, as defined 
above. Then, for X > fye 2 + An^/2e, the sequence of approximations defined 

in section 4 converge exponentially to f , i.e., 


max \f{x)~ f (M ’ XM] (x)\ <K x g M + ^jC, M , where C = ^9 0 /eX < 1. 

Here K x and K 2 are constants, independent of M , g is related to the lengths of the 
intervals and the singularities of /, (see equation (51)) , ( is relate d to the 

jumps in the derivatives of f (see equation (54)), and 6 0 = 1/ (2^1 - 1/(4A 2 )) . 

Before we illustrate the results of this theorem, we make two observations. First, if 
the jumps satisfy a bound of the form of (54) and fir < 1, then it follows that 

THEOREM 3 holds with <4 — Bo/eX. (This conclusion follows from the estimates derived 
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above and the fact that (7r£) M+1 — » 0, as M — » oo.) Secondly, we observe that if the 
jumps in satisfy a bound of the form | j/W(x 5 )j | < Af k , i.e., of the form of equation 
(54), but without the factor A:! present, then, again, Theorem 3 holds with ( = 0 o feX. 
(This follows from the fact that, for this case, the summation” in equation (53) is 
bounded by a constant times expfy^), as m — M — * oo.) 

We now summarize these observations in the following corollary. 

Corollary 4. If the jumps f^(x s ) J satisfy a bound of the form of (54) and 

fir < 1, or if they satisfy a bound of the form |[/ (fc) (x s )j| < Af k , s = 1,2, ...n, where 
A and f are positive constants, independent of k and s, then the results of Theorem 
3 hold with f = 9 0 jeX. In either of these cases, C be strictly less than one if 
X > ( l / 2 e ) y / l + e 2 = 0.533. 


6. More examples. In this section we consider three examples which serve to 
illustrate some of the convergence properties discussed in the previous section. 

Example 3: Let w be a non-integer constant and define f(x ) = cos(u;(:r — 7r)) , 
for — 7 r < x < 7 r, and f(x + 2tt) = f(x) for all x. We note that, since the cosine function 
is analytic in the entire finite part of the complex plane, the only singularities in / arise 
from the fact that cos (a; (a: — tt)) is not 2 tt— periodic. In this case, we set n — 1, with 
x 0 = — tt and X\ — i r. Then the quantity f] in equation (52) can be made arbitrarily small 
and, hence, the convergence of our approximations is dominated by the second term on 
the right side of equation (57). Also, for k = 0, 1,2, ..., / w has only one singularity in 
the interval (— tt, 7t] and this singularity lies at x = 7r, with the corresponding jumps 
given by 


(58) 



(— 1 ) 1+fc / 2 u; fc (1 — cos(2u/tt)) , 
( — l)( fc — 1 )/ 2 oa* sin(2o;7r), 


if k is even, 
if k is odd. 


Thus, the jumps in satisfy a bound of the form specified in equation (54) with 
£ = us, but without the factor of k\. Thus, if follows from COROLLARY 4 that our 
approximations converge exponentially to / at a rate which is bounded by C — fy/ eX. 
The Fourier coefficients of / are given by 

(-l) J+1 a;sin(2c<;7r) t (-l) J+1 j(l — cos(2w7r)) 

(59) = ir(j 2 - u 2 ) ’ b] = n(f - c 2 ) 

while the coefficients and bf I] are defined by 

(o) _ ( 2 fc+i) _ (2k+2) _ (~i) J ~ l ~ 1 k ? sin(2o;7r) ^-j. f 4a; 2 — (2 i + 1) 2 A 

a } -a,, a 3 - - n ^2_ ljJ 2) U ^ Aj 2 - (2i + l) 2 J \ 


(60) 


,(2 k) _ ,(2fc+i) (~l) j+1 i(l -cos(2a>7r)) k 

j ~ j n(j 2 -u 2 ) fi 


4a; 2 — (2 i + 1) 2 \ 

4i 2 _ (2 i + 1)2 ) > 
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for k = 0, 1,2, ... . Following the same arguments given in section 5, we find for this 
example, 




< 


K 

M 



R ( , 1 
M \2eJ\ 2 - 1/4 


M 


In Figure 9 we have plotted the normalized errors M (2e^\ 2 — 1/4) E <M ' XM \x) for 
A = 1 and several values of M, for the representative case when u = 7/5. In Figure 10 

we have plotted the normalized L^— errors M (2e yj\ 2 — 1/4) max |£’( Af,AM )(x)| as a 
function of M for several values of A, also for u> = 7/5. We shall discuss these results in 
section 7. 

Example 4: We let e be a positive constant and define f(x) = 1 j[x 2 + e 2 ), for 
— it < x < i r, and f(x + 2ir) = f{x), for all x. For this example, / has singularities in 
the complex plane at x = He, which he “close to” the interval [— 7r, 7 r], if e is small. 
Following the remarks before Theorem 3, for this example we set n = 2 and define 
xq — —it, xi = 0, and x 2 = n. Then we can select a\ = —n and 0.2 = 7r and find 

7T 

^ \At 2 + e 2 

In addition, each f^{x) has a singularity in the interval ( — 7r, 7 r], which lies at x = X 2 = 
7 r, with the corresponding jumps given by 



/ W (< = 


2k\ 




E (-i)'i 


r k+i x 

2j+h 


W k ~ 2 ^ 2 i 


(62) 

and 
it fol 


(7 r 2 + e 2 ) fc+1 ^ 

2A:! 

^—^jsin^ + litan- 1 ^)), * odd, 


/“'Wl = 0, if k is even. Using the fact that Jsin(^r) | < 1 for all real values of 2, 
ows from equation (66) that 


(63) 


/ W M' 


< 


tyfn 2 + e 2 \ a/tt 2 + c 2 


fc!, 


for all values of k. Thus, the jumps in satisfy a bound of the form specified in 
equation (54), with A = 2/(e\/7r 2 + e 2 ) and £ = (7 r 2 + c 2 ) -1 ^ 2 . Since < 1 for any 
real value of e, it follows from Corollary 4 that our approximations again converge 
exponentially to / at a rate bounded by the right side of equation (57), with 77 given by 
equation (61) and £ = 9 0 /e\. Thus, for “small” values of e, the rate of convergence is 
dominated by the value of 77, while for “large” values of e it is dominated by the value 

of C- 

To illustrate the rate of convergence of our approximations for this example, in 
Figure 11 we have plotted the quantities 


(64) 


Cm = (: 


max 


)(x)|)/( max | B (M-2,A(Jtf-2» (x) |)] 
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1/2 


as a function of 1/M for A = 1 and various values of c between 0.1 and 2. As M -► oo, 
approaches the asymptotic rate of convergence of our approxim ations. For the 
“small” value of e = 0.1, the asymptotic rate of convergence r) = n/ \Ar 2 + ^ — 0.999 
is indicated by the short dotted line near Cm = 1, while for the “large” value of e = 2, 
the asymptotic rate C = 0 o /e A = 0.212 is also indicated by a short dotted line. We note 
that, in both of these “extreme” cases, the figure seems to indicate that the estimate of 
the rate of convergence predicted by the theory is overly conservative, with the actual 
rate of convergence being smaller (better) in each case than the theory predicts. We 
shall discuss this example further in the next section. 

Example 5: Finally, let e be a real, positive number and define the function 
j-f x ) — ((tt 4- e ) 2 — x 2 ) -1 , for —n < x < 7r, and f(x + 2tt) = /(x), for all x. For this 
example, / has singularities at x = ±(tt + e), which lie “close to” the interval [— vr, tt], 
if e is small. For this example, we set n = 1 and define x 0 = -tt and x x = tt, and also 
set q:i = 0. Then 

( 65 ) * = VVe 


In addition, each f (k) (x) has a singularity in the interval (-tt, n], which lies at x x x 
7r, with the corresponding jumps given by 


/ (fe) w] 


( 66 ) 


2k\ (fc ^l /2 fk + V 

(e(2i t + e)) fc+1 j = o \2 j + 1, 

k\ { 1 1 


7 r 


k-2j 


(n + e) 


2 3 


7T 


+ € 1 e k+l (2tt + e) fc+1 


, k odd, 


and /^ fc) ( 7r )] = 0, if k is even. Using equation (66), we find 


(67) 




and, hence, the bound (54) holds with C = Thus, for this example, si nce ftr > 1 
when € < 7T, Theorem 3 guarantees exponential convergence only for A > \/e 2 + 47 r £ /2e 
7r/ee = 1. 1557/e, as e — ► 0. Alternatively, fo r a fixed value of A, the requirement that 
£ < 1 will be satisfied only if e > e 0 = ^/(e^A 2 - 1/4), with e 0 = 1-335 when A = 1. 

To illustrate the rate of convergence of our approximations for this example, in 
Figure 12 we have plotted the quantities (m defined by equation (64) as a function 
of 1/M for A = 1 and e = 0.25, 0.375, 0.5, and 1. As M -* oo, Cm approaches the 
asymptotic rate of convergence of our approximations. It is obvious from the plot that 
the lower bound on e obtained above (e 0 = 1.335) is too conservative, with the actual 
value of e 0 being closer to 0.4 (for A = 1). We shall discuss these results further in the 

next section. 

7. Discussion and Conclusions. We now make some observations concerning 
the convergence properties of the approximations we have constructed, particularly m 
light of the examples we have considered. 
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For the function of Example 1, each Fourier coefficient is just a monomial in (i/i). 
i.e., bj = -2/j. In fact, this function is just a constant multiple of the function h 0 , 
which was the basic “building block” for the convergence proofs of section 5. In this 
case, explicit expressions for the Fourier coefficients b ^ (equation (30)) were obtained, 
where their 0(1/ j M + 2 ) behavior is evident. Since /<*>(*) = 0, for k> 2, the “need” to 
include the basis functions S k with k > 2 in the expressions for R k is due solely to the 
jump in the k th derivative of S t , for l < k (see, e.g., equation (23)). This observation 
suggests that it might be worthwhile to try to construct a new set of basis functions, say, 

| iSw , which have the property that S k has a jump in only its k derivative at x 0. 

Such functions can be constructed in a straightforward maimer, by expressing S k as an 
infinite linear combination of the S h with l >k. However, this series representation of 
S k fails to converge for all x in the range -7r < x < n, and hence is not useful. Even if 
these expressions for S k did converge over the entire interval, it is not clear that they 
would have any practical advantage over the f mictions {S k } which we have introduced. 

The function of Example 2 is a piece-wise polynomial function, as defined in Lemma 
1, with m = 2, and has n = 4 points of singularity in the interval -n < x < n. 
Consequently, each of its Fourier coefficients is a polynomial in (1/j). Even though 
this function has multiple singularities, the approximations we constructed are able to 

handle them in a straightforward manner. 

Each Fourier coefficient of the function of Example 3 can be expressed as an infinite 
series in powers of (1/j). However, the construction of our sequence of exponentially 
convergent approximations follows exactly the procedure as in the polynomial case. In 
fact, the only technical difference is the inclusion of the jump in f (k) in the definition 
of the constants A kJ (see equation (21)). For this example, the function / is an entire 
function, with no singularities in the finite part of the complex x— plane. In particular, 
it is analytic on [— tt, tt], but is not 27r-periodic. Consequently, the only “singularities” 
in the 27r-periodic extension of / occur at x = tt. Because of these properties, the 
rate of convergence of our approximations for this example is the same as the rate of 

convergence for the function ho, i.e., f — 9o/eX. 

Each Fourier coefficient of the function / of Example 4 is again an infinite series in 
(1/j), although these coefficients must now be computed numerically. As in Example 
3, / is analytic at each value of x in the interval [— tt, 7t], but it is not 27r-periodic. 
However, for this case, / has (purely imaginary) singularities in the finite part of the 
complex plane at x = ±ie. As a consequence, the rate of convergence of the approxima- 
tions depends upon the value of e, with the rate of convergence approaching 1 as e -► 0. 
However, as e —* oo, the rate of convergence is again the same as that for the function 
ho- Also, for any real, positive value of e, the approximations converge exponentially for 

any value of A > Ao = (l/2e)y/l + e 2 - 

The function of Example 5 is also analytic on [— 7r, 7 t], but, again, is not 27r-periodic. 
In this case, however, f has real singularities at x = ±{tt +e), and the rate of convergence 
of our approximations again depends upon e. In this case, however, for sufficiently small 
values of e the value of A must be greater than a specific lower bound, which depends 
upon e and which increases as e -» 0, in order to ensure exponential convergence. 
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In comparing the results of Examples 4 and 5, we note that, in both cases, there are 
singularities of / which lie a distance e from the interval [ — tt, 7 t]. In Example 5, when 
e is “small”, these singularities lie close to the points where the 27r-periodic extension 
of / has real singularities, namely, at x = ±7r. In this case, the presence of these “near- 
by” singularities seems to adversely affect the convergence of our approximations, in 
the sense that we are forced to require that A increase as e decreases. By contrast, in 
Example 4 the singularities lie “far” from the points where the 27r-periodic extension 
of / has real singularities, even when e is small. In this case, although the rate of 
convergence is affected by the value of e through the parameter 77, it is not necessary to 
require that A increase as e becomes smaller. Thus, the “nearness” of the singularities 
of / to the points {x.,} seems to be of more significance than simply the nearness of the 
singularities to the interval [— n, 7r] . 

In the formulas above, we have consistently used (because of its simplicity) the 
simple (loose) bound (9 0 in our estimates of the rates of convergence of our approxi- 
mations. These estimates could be improved somewhat by using instead the sharper 
estimate derived in the Appendix. For example, the plots in Figure 8 of the Loo-errors 
for Example 1, which were normalized using 60 , clearly show, by the downward trends 
as M — » co, that the actual rate of convergence is smaller than the rate d 0 /\e which 
was used in the normalization. In Figure 13, we have plotted the maximum (L 0 c ) 
errors max |E'( Af,AM )(x)| for Example 1, normalized by M(Xe/D 0 9o) M , which involves 
the sharper error bound derived in the Appendix, where Do is defined by equation 
(73). Each of these plots are clearly bounded as M becomes large, and appear to be 
approximately linear as M — > 00, with a slope that approaches zero as A becomes 
large. These observations are consistent with an asymptotic rate of convergence given 
by Z>o#o/Ae < 9o/Xe. 

As mentioned in the introduction, Gottlieb and several of his co-investigators 
[3], [4], [5], [6], [7] have proposed and analyzed a technique using the Gegenbauer polyno- 
mials which also leads to a series which converges exponentially to / in the maximum 
norm (and, hence, also eliminates the Gibbs phenomena). Their technique requires a 
knowledge of only the first 2 N + 1 Fourier coefficients of /, whereas our method requires 
a knowledge of the location and magnitude of the jumps of / and its first N derivatives 
within the interval [-7T,7rj. Consequently, since our technique requires more “informa- 
tion” about the function than Gottlieb’s method, it is not surprising that our approxima- 
tions converge much more rapidly than his do. In particular, the error in his approxima- 
tions contains a term proportional to N 2 (qr) N , where q T = (fi + 2a)P(27rea)~ a P~ 0 < 1, 
and a and ft can be chosen so that the convergence rate is optimized. The correspond- 
ing term in our approximations (see THEOREM 2) is N~ 1 (9q/ e) N , where we have set 
A = 1 and, hence, M — N. For the function / of our Example 1 (which is also the 
first example in [3]), these terms are (approximately) 7V 2 (0.4763) W and ./V _1 (0.2124) N , 
respectively. In Figure 14 we have plotted log 10 max |£A w,w )(x) as a function of N 
for Example 1 and Example 3 using our approximations and Gottlieb’s approximations. 
The difference in the rates of convergence is apparent from the figure 

We should also point out that in [3] the assumption was made that p > 1, where 
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p is the distance of the nearest singularity of / from the basic interval [—1,1]. (The 
corresponding assumption for our analysis would be that p > 7r.) This condition is 
relaxed somewhat in [5]. However, we impose no such restriction on our functions /. 

Several other methods have been proposed to accelerate the convergence of Fourier 
series and/or to eliminate (or at least mitigate) the Gibbs’ phenomena. These methods 
include the construction of Fourier-Pade approximations [2] (mentioned in the intro- 
duction), Fejer’s method [1], based on Cesaro sums, and Lanczos’ “sigma” method [11]. 
In Figure 15, we have compared these three approximations for the function f(x) = x 
of Example 1, along with the original Fourier series partial sum F (4) (equation (4)), 
Gottlieb’s approximation [3], and our approximation /^ ,4 ^ (equation (14)). It is inter- 
esting to note that, although the approximations of Fejer and Lanczos have essentially 
eliminated the Gibbs phenomena, these approximations are of poor quality in an L 2 
sense. The simple reason for this is that each of these methods retains the same form of 
the approximation and just replaces the original Fourier coefficients with certain other, 
modified coefficients. It is well known (see [1], for example) that, for any linear combi- 
nation of a finite number of sine and cosine functions, the original Fourier coefficients 
are optimal in an L 2 sense. Thus, the approximations of Fejer and Lanczos can be no 
better in the L 2 sense than the original Fourier series partial sum. 

Currently, we are investigating several topics related to the approximations we have 
proposed, and they will be discussed in future papers. These topics include the problem 
of determining estimates of the locations and the magnitude of the jumps in / and its 
derivatives in the interval [— 7T, 7 r] directly from the Fourier coefficients of /. Such esti- 
mates will allow the basis functions defined here to be used in a manner similar to that 
described above, when only a finite number of the Fourier coefficients of / are known. 
Some preli mi nary investigations have indicated that the new approximations will still 
retain the property of exponential convergence, although the rate of convergence of the 
new approximations will undoubtedly not be as rapid as the approximations described 
here. Applications of these ideas directly to the solution of differential equations (par- 
ticularly, boundary value problems) are also being investigated, as well as an extension 
of the ideas presented here to other orthogonal basis (eigenfunction) expansions. 
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Appendix: Derivation of 6q 

In this appendix, we first present a brief derivation of the bound 


( 68 ) 


* i = 0 


1 - 


2i+ 1 
k 2 j 


< K6 ", 0 O (A) = 


2 v /l-l/(4A2)’ 


where M = 2m + 1 and j > AM + 1, which was used in connection with the discussion 
of the function of Example 1, as well as the function /to which was used in the proof of 
Theorems 2-3. To show that equation (68) holds, we first observe that 


2i + 1 ^ 2m + 1 ^ J_ 

2j ~ 2(AM + 1) 2A’ 

for 0 < i < m and j > AM + 1. Using this result, we can write 



Using this bound in equation (68) we find, for all A > A > 1/2, 

(69) o M <ke^, K = (\/l - 1/(4A 2 )^ . 


In the special case when AM is an integer, we can obtain a sharper bound for 6 M . 
(Recall that M is an odd integer.) To do this, we first define J = 2(AM + 1) and then 
note that 



-.o \ u / 

M+l (J-M- 2)!! 
(J + M)!! ' 


(70) 
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(Here we have used that usual notation k\\ = 1 - 3 ■ 5 - • • fc, for any odd integer k.) 
Then, using Stirling’s approximation to the factorial, we find that k\\ — ► y/2 k(k/e) k/2 , 
as k —*■ oo. Using this approximation in equation (70), we find that 

(71) e M < K • ( D 0 9o) M , for A > A > 1/2, 


where 

(72) 

and 

(73) 


k = k{ a) = 


'2A + 1 


_J ] g 1 /* . < — 

2X-1J yi - 1/(4A 2 ) \/l — 1/(4A 2 ) 


Dq = -Do(A) = e 


2A — 1 
2A + 1 


1, as A — » oo. 


We remark that, although our derivation is valid only for the case when AM is 
an integer, the final result (equations (71)-(73)) holds for any positive value of A. The 
reason for this is that the expression for k\\ can be written as (fc!)/(2 (fc_1)/2 ((A: - l)/2)!), 
and this expression can be defined for noninteger values of k in terms of the Gamma 
function. Since the estimate we used for k\\ (following equation (70)) also holds for 
noninteger values of k, it follows that the expressions (71)- (73) hold for noninteger 
values of AM as well. 
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2) The function / (solid line), the function R x (dashed line), and the difference 
fi = / — Ri (dotted line), for Example 1. 
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3) The function / (solid line), the function (the barely visible dashed line), 
and the magnified error 50(/ — /^ ,3 ^) (dotted line), for Example 1. 
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4) The function / for Example 2 (solid lines), along with the Fourier series partial 
sum of / (see equation (4)) with N = 10 (dashed line), and the error / — F 
(dotted line). 
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6) The function / (solid line), the function (the barely visible dashed line), 

and the magnified error 50(/ - f( l &) (dotted line), for Example 2 
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7) The normalized errors M(\e/6 0 ) M E^ M,XM \x) for Example 1, when A = 1 and 
M = 5 (solid line), M = 9 (dashed line), and M = 13 (dotted line). Note the more 
or less “uniform” distribution of the error over the interval [— 7r, 7r] (especially as M 
increases), as opposed to the error in the original Fourier series partial sums, which is 
highly concentrated near x = 0. 









9) The normalized errors M(Xe/9 0 ) M £ (M ’ AM) (*) for Example 3, when A = 1 and 
M = 5 (solid line), M = 9 (dashed line), and M = 13 (dotted line). Note the more 
or less “uniform” distribution of the error over the interval [ tt, 7t] (especially as M 
increases), as opposed to the error in the original Fourier series partial sums, which is 
highly concentrated near x = tt. 







11) The quantities (m = (max j J E7 (JV/ ’ AA ' / > (^) |) / (max E^ M ~ 2, ^ M ~ 2 ^{x) for 
Example 4 plotted as a function of 1/M for A = 1 and for e = 0.1, 0.25, 0.375, 0.5, 1. and 
2. As M — > oo, Cm approaches the asymptotic rate of convergence of our approximations. 
The theoretical asymptotic rates of convergence 77 = 7 r/\Ar 2 + e 2 = 0.999 (for e = 0.1) 
and C = #o/ e A = 0.212 are each indicated by a short dotted line. 




12) The quantities Cm = [(max |e (m ' am) (x)|) / (max | E^ M 2>A(M 2)) ( x )|)] fol 
Example 5 plotted as a function of 1/M for A = 1 and for e = 0.25, 0.375, 0.5, and 1 
As M -MX), Cm approaches the asymptotic rate of convergence of our approximations 
The critical value Cm = 1 is indicated by the short dotted line. 













l°g n n [max I error | ] 



14) A plot of the log 10 {max |Errorj) for the functions of Example 1 (dotted lines) 
and Example 3 (solid lines) using Gottlieb’s approximation [3] (diamonds) and the 
approximations equation (22) (circles), plotted as a function of N. 
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15) A comparison of the function f(x) = x of Example 1 (denoted by circles) 
and the original Fourier series partial sum (equation (4)), the Fourier-Pade ap- 
proximation F 3 i [2], Fejer’s approximation [1], Lanczos’ “sigma” approximation [11], 
Gottlieb’s approximation [3], and the approximation (equation (14)). All of these 
approximations are based on only the first 4 terms in the original Fourier partial sum. 
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